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ABSTRACT 

Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how 
the measure of star formation is affected by the spatial scale at which a galaxy is observed. 

Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and 
other directly associated physical parameters such as the attenuation. 

Methods. We have carried out a multi-scale, pixel-by-pixel study of the nearby galaxy M33. Assembling FUV, Her, 8 jum, 24 jum, 
70 jum, and 100 jum maps, we have systematically compared the emission in individual bands with various SFR estimators from a 
resolution of 33 pc to 2084 pc. 

Results. We have found that there are strong, scale-dependent, discrepancies up to a factor 3 between monochromatic SFR estimators 
and lla+24 pm. The scaling factors between individual IR bands and the SFR show a strong dependence on the spatial scale and on 
the intensity of star formation. Finally, strong variations of the differential reddening between the nebular emission and the stellar 
continuum are seen, depending on the specific SFR (sSFR) and on the resolution. At the finest spatial scales, there is little differential 
reddening at high sSFR. The differential reddening increases with decreasing sSFR. At the coarsest spatial scales the differential 
reddening is compatible with the canonical value found for starburst galaxies. 

Conclusions. Our results confirm that monochromatic estimators of the SFR are unreliable at scales smaller than 1 kpc. Furthermore, 
the extension of local calibrations to high redshift galaxies presents non-trivial challenges as the properties of these systems may be 
poorly known. 

Key words. Galaxies: individual: M33; galaxies: ISM; galaxies: star formation 


1. Introduction 

As we observe galaxies across the Universe, their evolution from 
highly disturbed proto-galaxies at high redshift to the highly or¬ 
ganised systems common in the zoo of objects we see in the 
nearby Universe is striking. One of the most important pro¬ 
cesses that drives this evolution is the transformation of the pri¬ 
mordial gas reservoir into stars, which form heavy elements that 
are ejected into the intergalactic medium during intense episodes 
of feedback. In other words, if we want to understand galaxy 
formation and evolution across cosmic times, we need to under¬ 


stand the process of star formation in galaxies. To do so, it is 
paramount to be able to measure star formation as accurately as 
possible. 

The most direct way to trace star formation is through the 
photospheric emission of massive stars with lifetimes of up to 
~ 100 Myr, which dominate the ultraviolet (UV) energy bud¬ 
get of star-forming galaxies. An indirect star formation tracer 
is the Ha recombination line (or any other hydrogen recom¬ 
bination line) from gas ionised by the most massive stars that 
are around for up to ~ 10 Myr. However, both the UV emis¬ 
sion and the Ha line are severely affected by the presence of 
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dust, which absorbs energetic photons and reemits their energy 
at longer wavelengths. From the inception of the far-infrared 
era with the launch of IRAS (Infrared Astronomical Satellite, 
Neugebauer et al. 1984), the emission of the dust has been used 
as a powerful tracer of star formation from local galaxies up to 
high-redshift objects, resulting in a tremendous progress of our 
understanding of galaxy evolution in general and of the physical 
processes of star formation in particular. 

The launch of the Herschel Space Observatory (Pilbratt et al. 

2010) has opened new avenues for the investigation of star for¬ 
mation in the far-infrared not only in entire galaxies, but also 
within nearby galaxies at scales where physical processes such 
as heating and cooling are localised. Herschel matches the angu¬ 
lar resolution of 5-6"of Spitzer in the mid-infrared (Fazio et al. 
2004; Rieke et al. 2004), and of GALEX in the UV (Galaxy Evo¬ 
lution Explorer, Martin et al. 2005). 

Yet, measuring local star formation in galaxies remains an 
important challenge. Eor instance, Kennicutt et al. (2007); Bigiel 
et al. (2008) found seemingly incompatible star formation laws 
with the same dataset. Such a difference could be due to the 
distinct ways star formation is measured in galaxies (Liu et al. 

2011) , different SER estimates leading to variations of 10-50% 
of the molecular gas depletion timescale (Leroy et al. 2012, 
2013). 

The measurement of star formation relies upon three main 
assumptions. 

- Eirst of all, a well-defined and fully sampled initial mass 
function (IME) is assumed. This is necessary to relate the 
measured power output from massive, short-lived stars to 
the total mass of the stellar population of the same age. Mas¬ 
sive stars only account for a minor fraction of the total mass 
of stellar populations, even in the youngest star-forming re¬ 
gions, which contain the highest proportion of such stars. 

- Star-formation-tracing bands need to be sensitive mainly to 
the most recent episode of star formation. Contamination 
from emission unrelated to recent star formation, such as ac¬ 
tive nuclei and older stellar populations, needs to be negligi¬ 
ble. 

- A well-defined star formation history. Too few star-forming 
regions would induce rapid variations of the SER with time. 

These assumptions, which are not exhaustive, may already be 
problematic for some entire galaxies (Boselli et al. 2009). At 
small scales, they are unlikely to hold true across an entire spiral 
disk. 

If we want to understand star formation laws in the era of 
resolved observations, it is therefore crucial to understand when, 
how, and from which spatial scale we can measure star formation 
reliably. In particular, we need to understand how star forma¬ 
tion tracers relate to each other in galaxies from the finest spatial 
scale, at which H ii regions are resolved, to large fractions of a 
spiral disk. Recent results show a systematic variation of star 
formation tracers with spatial scale, which could be due to the 
presence of diffuse emission unrelated to recent star formation 
(Li et al. 2013): ~ 20 - 30% of the EUV luminosity from a 
galaxy is due to stars older than 100 Myr (Johnson et al. 2013; 
Boquien et al. 2014) and a 30% to 50% of Ha is diffuse (Thilker 
et al. 2005; Crocker et al. 2013). Measuring local star forma¬ 
tion is made even more difficult by the fact that indirect tracers 
of star formation (the ionised gas and dust emission) may not be 
spatially coincident with the direct tracer of star formation, the 
UV emission (Calzetti et al. 2005; Relano & Kennicutt 2009; 
Verley et al. 2010; Louie et al. 2013; Relano et al. 2013). Such 
offsets can also be seen in the Milky Way in NGC 3603, Carina, 


or the OB associations in Orion for instance. This questions the 
real meaning of SER measurements at local scales. 

These offsets along with other processes such as stochastic 
sampling of the IME or the insufficient number of star-forming 
regions and/or molecular clouds in a given region could be one of 
several reasons for the observational breakdown of the Schmidt- 
Kennicutt law on scales of the order of -100-300 pc (Calzetti 
et al. 2012), which has been found in M33 by Onodera et al. 
(2010) and Schruba et al. (2010). The complex interplay be¬ 
tween various processes at the origin of the breakdown of the 
Schmidt-Kennicutt law at small spatial scales has recently been 
analysed by Kruijssen & Longmore (2014). 

With the availability of resolved observations of high- 
redshift objects with ALMA and the JWST by the end of the 
decade, understanding whether and how we can measure star 
formation at local scales is also of increasing importance. We 
address this question through a detailed study of star formation 
tracers at all scales in the nearby late-type galaxy M3 3. Thanks 
to its proximity (840 kpc, corresponding to 4.07 pc/", Ereed- 
man et al. 1991), relatively low inclination (56°, Regan & Vo¬ 
gel 1994), and large angular size (over 1° across), M33 is an 
outstanding galaxy for such a study. It has been a popular tar¬ 
get for a large number of multi-wavelengths observations and 
surveys in star-formation tracing bands from the EUV with the 
GALEX Nearby Galaxies Survey (NGS, Gil de Paz et al. 2007), 
to the EIR with Herschel in the context of the HerM33es survey 
(Kramer et al. 2010), including Spitzer mid-infrared data (Ver¬ 
ley et al. 2007) as well as Ha narrow-band imaging (Hoopes & 
Walterbos 2000). 

In Sect. 2 we present the data, including new observations 
recently obtained by our team, and how data processing was car¬ 
ried out. We compare various SER estimators at different scales 
in Sect. 3. We examine in detail the properties of dust emission 
with scale in order to measure the SER from monochromatic in¬ 
frared bands in Sect. 4. We investigate the relative fraction of 
attenuated and unattenuated star formation with scale in Sect. 5. 
Einally, we discuss our results in Sect. 6 and we conclude in 
Sect. 7. 

2. Observations and data processing 

2.1. Observations 

To carry out this study, we consider the main star formation trac¬ 
ers used in the literature: the emission from young, massive stars 
in the EUV, the ionised gas recombination line Ha, and the emis¬ 
sion of the dust at 8 pm, 24 pm, 70 pm, and 100 pm. We do not 
explore dust emission beyond 100 pm as prior work has shown 
that longer wavelengths are not good tracers of star formation 
(Bendo et al. 2010, 2012; Boquien et al. 2011) and the scale 
sampled with Herschel becomes coarser. We also forego radio 
tracers as they are not as widely used. 

The EUV GALEX data from NGS were obtained directly 
from the GALEX website through galexview^ The observation 
was carried out on 25 November 2003 for a total exposure time 
of 3334 s. 

HQf-i-[Nii] observations were carried out in November 1995 
on the Burrel Schmidt telescope at Kitt Peak National Observa¬ 
tory. They consisted in 20 exposures of 900 s, each covering a 
final area of 1.75 x 1.75 deg^. This map has been continuum- 
subtracted by scaling an off-band image using foreground stars. 
The observations and the data processing are analysed in detail 
in Hoopes & Walterbos (2000). 

^ http://galex.stsci.edu/GalexView/ 
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The Spitzer IRAC 8 pm image sensitive to the emission of 
Polycyclic Aromatic Hydrocarbons (PAH) and the MIPS 24 pm 
image sensitive to the emission of Very Small Grains (VSG) 
were obtained from the NASA Extragalactic Database and have 
been analysed by Hinz et al. (2004) and Verley et al. (2007). 

The PACS data at 70 pm and 100 pm, which are sensitive to 
the warm dust heated by massive stars, come from two different 
programmes. The 100 pm image was obtained in the context 
of the Herschel HerM33es open time key project (Kramer et al. 
2010, observation ID 1342189079 and 1342189080). The ob¬ 
servation was carried out in parallel mode on 7 January 2010 
for a duration of 6.3 h. It consisted in 2 orthogonal scans 
at a speed of 20'7s, with a leg length of 70'. The 70 pm 
image was obtained as a follow-up open time cycle 2 pro¬ 
gramme (OT2_mboquien_4, observation ID 1342247408 and 
1342247409). M33 was scanned on 25 June 2012 at a speed 
of 20"/s in 2 orthogonal directions over 50' with 5 repetitions 
of this scheme in order to match the depth of the 100 pm im¬ 
age. The total duration of the observation was 9.9 h. Reduced 
maps are available on the Herschel user provided data product 
website^. 

2.2. Additional data processing 

The GALEX data we obtained from galexview were already 
fully processed and calibrated, we therefore did not carry out 
any additional processing. 

We corrected the continuum-subtracted Ha map for [Nii] 
contamination, which according to Hoopes & Walterbos (2000) 
accounts for 5% of the Ha flux in the narrow-band Alter. We 
have also removed subtraction artefacts caused by bright fore¬ 
ground stars. To do so we have used irae’s imedit procedure, 
replacing these artefacts with data similar to that of the neigh¬ 
bouring background. 

The Spitzer IRAC and MIPS data we used were processed in 
the context of the Local Volume Legacy survey (EVE, Dale et al. 
2009). No further processing was performed. 

Even though in the context of the HerM33es project we al¬ 
ready reduced and published 100 pm data (Boquien et al. 2010b, 
2011), these observations were processed with older versions 
of the data reduction pipeline. To work on a fully consistent 
set of Herschel PACS data and to take advantage of the recent 
improvements of the pipeline, we have reprocessed the 100 pm 
from the HerM33es survey along with the new 70 pm data. To 
do so we have taken the raw data to level 1 with HIPE version 9 
(Ott 2010), flagging bad pixels, masking saturated pixels, adding 
pointing information, and calibrating each frame. In a second 
step, to remove the intrinsic 1// noise of the bolometers and 
make the maps, we used the Scanamorphos software (Roussel 
2013), version 19. We present the new 70 pm map obtained for 
this project in Eig. 1. 

2.3. Correction for the Galactic foreground extinction 

To correct the EUV and Ha fluxes for the Galactic foreground 
extinction, we have used the Cardehi et al. (1989) extinction 
curve, including the O’Donnell (1994) update. We have as¬ 
sumed E(B - V) = 0.0413 as indicated by NASA/IPAC Infrared 
Science Archive’s dust extinction tool from the Schlegel et al. 
(1998) extinction maps. This yields a correction of 0.34 mag in 
EUV and 0.11 mag in Ha. 

^ http://WWW.cosmos.esa.int/web/herschel/ 
user-provided-data-products 


2.4. Astrometry 

To carry out a pixel-by-pixel analysis, it is important that the rel¬ 
ative astrometric accuracy of ah the bands is signiflcantly better 
than the pixel size. A first visual inspection reveals a clear offset 
between the new 70 pm data we present in this paper and the 
100 pm data presented in Boquien et al. (2010b, 2011). When 
comparing the 70 pm and 100 pm images with the Ha image of 
Hoopes & Walterbos (2000), we found that the 70 pm map corre¬ 
sponded more closely to the Ha emission across the galaxy and 
was consistent with data at other wavelengths. We have there¬ 
fore decided to shift the 100 pm image to match the 70 pm map 
astrometry. To determine the offset, we have compared the rela¬ 
tive astrometry of the 160 pm images obtained in the context of 
HerM33es and the OT2_mboquien_4 programme. As the 160 pm 
is observed in parallel with the 70 pm or the 100 pm, they have 
the same astrometry. The offset between the 160 pm maps be¬ 
tween these two programmes is therefore the same as the offset 
between the 70 pm and the 100 pm maps. We have determined 
an offset of ~ 5" (4.83" towards the East and 1.25" towards the 
North) and applied this to the 100 pm image. When comparing 
the corrected 100 pm band with the 8 pm and 24 pm images 
we can see small region-dependent offsets of the order of 1-2". 
The variation of this offset from one region to another leads us 
to think that at least part of it reflects physical variations in the 
emission of the various dust components in M33. In addition, as 
we will see later, we carry out this study at a minimum pixel size 
8", much larger than any possible systematic offset. We con¬ 
clude that the relative astrometry of our images is sufhcient to 
reach our goals. 

2.5. Pixel-by-pixel matching 

Because pixel-by-pixel analysis will be central for this study, it 
is crucial to match ah the images to a common reference frame. 
To do so, it is important that ah bands share a common point 
spread function (PSE). To ensure this, in a first step we have 
convolved all the images to the PACS 100 pm PSP using the 
dedicated kernels provided by Aniano et al. (2011). We have 
then registered these images to a common reference frame with a 
pixel size ranging from 8", slightly larger than the PACS 100 pm 
PSP, to 512", by increments of 1" in terms of pixel size, using 
irae’s wregister procedure with the drizzle interpolant. This al¬ 
lows us to sample ah scales from fractions of H ii regions at 33 pc 
(8") to large fractions of the disk at 2084 pc (512"). The upper 
bound is limited by the size of the galaxy. Increasing to larger 
physical scales would leave us with too few pixels in M3 3. We 
present some of the Anal, convolved, registered, and background 
subtracted maps used in this study for a broad range of pixel 
sizes in Pig. 2. 

To compute flux uncertainties, we have relied on the 33 pc 
scale maps. We took into account the uncertainty on the back¬ 
ground determination, which is due to large scale variations, and 
the pixel-to-pixel noise. The former was measured as the stan¬ 
dard deviation of the background level measured within 10 x 10 
pixel square apertures around the galaxy using irae’s imexamine 
procedure. The latter was measured as the mean of the standard 
deviation of pixel fluxes in these apertures around the galaxy. 
We then summed these uncertainties in quadrature. Por maps at 
lower resolution, we simply scaled the uncertainties on the back¬ 
ground with the square of the pixel size, and the pixel-to-pixel 
uncertainties with the pixel size. Direct measurements on lower 
resolution maps yielded uncertainties consistent with the scaled 
ones. 


Article number, page 3 of 15 



+ 31°00'00.0" 


.0.0150 

0.0135 

0.0120 

0.0105 

0.0090 

0.0075 

0.0060 


50'00.0" 


. • > 


o 

o 

o 

CM 

u 

0 

Q 


40'00.0" 




Jt Ti 


30'00.0" 


» 1 ^ 


0.0045 


- - 0.0030 


o 

0 

CA) 

0 

0 

0 


0.0015 


+ 30°20'00.0" 


30.00s 35m00.00s 30.00s 34m00.00s 30.00s 

RA a2000) 


33m00.00s 30.00s 


' 0.0000 
Ih32m00.06s 


Fig. 1. Map of M33 at 70 jim obtained with the Herschel PACS instrument in the context of a cycle 2 programme (OT2_mboquien_4, observation 
ID 1342247408 and 1342247409; the original map is available from the link given in footnote 2). The image is in Jy/arcsec^ and the colours follow 
an arcsinh scale indicated by the bar on the right side of the figure. The physical scale is indicated by the white line in the bottom-right corner of 
the figure, representing 1 kpc. Each pixel has a size of 1.4". 


2.6. Removal of the stellar pollution in infrared bands 

In a final data processing step, we have removed the stellar con¬ 
tamination in the 8 pm and 24 pm bands. To do so, we have 
assumed that the Spitzer/IRAC 3.6 pm image is dominated by 
stellar emission, following the analysis of Meidt et al. (2012). 
We have then scaled this image to predict the stellar emission at 
8 pm and 24 pm and subtracted it from these images. We have 
assumed a scaling factor of 0.232 at 8 pm and 0.032 at 24 pm, 
following Helou et al. (2004). We should note that this scaling 
factor can change quite significantly with the star formation his¬ 
tory (e.g. Calapa et al. 2014; Ciesla et al. 2014). 


3. Comparison of SFR estimators at different scaies 

3.1. Presentation of monochromatic and composite SFR 
estimators 

Ideally, a good SFR estimator has a solid physical basis and 
is devoid of biases. Thus, because they trace directly or indi¬ 


rectly the emission from young, massive stars, the attenuation^- 
corrected FUV or Her should in principle be ideal estimators. In 
practice however, the presence of biases is a real problem since 
it shows that other factors unrelated to recent star formation can 
contribute to the emission in star-formation tracing bands. For 
instance, in the case of monochromatic IR tracers, such factors 
are the contribution from old stars, changes in the opacity of the 
ISM (interstellar medium), or in the IR SED (spectral energy 
distribution). 

When no attenuation measurement is available, a popular 
method developed over the last few years has been to combine 
attenuated and attenuation free tracers (Calzetti et al. 2007; Ken- 
nicutt et al. 2007; Leroy et al. 2008; Kennicutt et al. 2009; Hao 
et al. 201 1). Unfortunately, how to combine such tracers remains 
uncertain. Calzetti et al. (2007) and Kennicutt et al. (2009) found 

^ We distinguish between the extinction which includes the absorption 
and the scattering out of the line of sight, and the attenuation which also 
includes the scattering into the line of sight. In practice in this study we 
only have access to the attenuation and not to the extinction. See for 
instance Sect. 1.4.1 of Calzetti (2013). 
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33 pc (S'Vpixel) 260 pc (64"/pixel) 2084 pc (512"/pixel) 



Fig. 2. Convolved images registered to a common reference frame 
at 33 pc (S'Vpixel, left), 260 pc (64"/pixel, centre), and 2084 pc 
(512"/pixel, right). Each row represents a different star formation trac¬ 
ing band, from top to bottom: GALEX EUV, Her, IRAC 8 yum, MIPS 
24 yum, PACS 70 yum, and PACS 100 jjim. Blue pixels have a low flux 
density whereas red pixels have a high flux density, following an arc- 
sinh scale. The colours used are simply chosen to best represent the 
large dynamical range of intensities across all bands and all pixel sizes 
and should be used in a qualitative sense only. 


different scaling factors when combining dust emission at 24 jim 
with Her, likely because of different scales probed: 500 pc for 
the former and entire galaxies for the latter, and therefore differ¬ 
ent timescales (Calzetti 2013). According to Leroy et al. (2012), 
the universality of composite estimators remains in doubt. One 
of the main issues comes from the diffuse emission and whether 
it is linked to recent star formation or not. In M3 3, the fraction 
of diffuse emission is high: 65% in FUV, and from 60% to 80% 
in the 8 yum and 24 yum bands, with clear variations across the 
disk for the latter two (Verley et al. 2009). While some methods 
have been suggested to remove the diffuse emission linked to 
old stars (Leroy et al. 2012), they rely on uncertain assumptions. 
We therefore cannot rely a priori on such tracers as an absolute 
reference. Yet, how monochromatic and hybrid SFR estimators 
compare may still yield useful information on star formation in 


Table 1. SFR estimators. 



Monochromatic 


Band 

log Cband k Method 

Reference 


FUV 

-36.355 

1.0000 

TheoreticaF 

1 

Ha 

-34.270 

1.0000 

TheoreticaF 

1 

24 pm 

-29.134 

0.8104 


2 

70 pm 

-29.274 

0.8117 

Ha'’ 

3 

100 pm 

-37.370 

1.0384 

Ha'’ 

3 


Hybrid 

Band 

log Cbandl 

kbandl-band2 

Method 

Reference 

Ha-\-24 pm 

-34.270 

0.031 

Rt? 

2 

FUV-r24 pm 

-36.355 

6.175 

Ha-\-24 pm 

4 


References. (1) Murphy et al. (2011), (2) Calzetti et al. (2007), (3) Li 
et al. (2013), (4) Leroy et al. (2008) 

Notes. Monochromatic: logESFR = log Cband + k x log Sband; Hybrid: 
log 2SFR — log Cbandl [^bandl kbandl-band2 ^ Sband2]5 with 2SFR 

in Mq yr"^ kpc“^, S defined as vSy in W kpc“^, and C in M© yr"^ W"k 
Empirical estimators have been calibrated on individual 
star-forming regions on typical scales of the order of 
-200-500 pc. 

^ Based on Starburst99 (Leitherer et al. 1999). ^ Extinction 
corrected, calibrated against near-infrared hydrogen 
recombinations lines (e.g., Patr or Bry). 

M3 3. We consider the restricted set of monochromatic and hy¬ 
brid SFR estimators presented in Table 1 . 

Before comparing these SFR estimators we have to add a 
word of caution. In some cases, especially at the smallest spatial 
scales, the concept of an SFR in itself may not be valid (for a de¬ 
scription of the reasons see Sect. 3.9 of Kennicutt & Evans 2012, 
in particular: IMF sampling, age effects, and the spatial exten¬ 
sion of the emission in star-formation-tracing bands in compari¬ 
son to the resolution). In the context of this study, IMF sampling 
is not likely to be a particular issue. A scale of 33 pc corresponds 
to the Stromgren radius of a 3000-5000 M©, 4-5 Myr old stellar 
cluster. Such a cluster would already be massive enough not to 
be too affected by stochastic sampling (Fouesneau et al. 2012). 
However, we cannot necessarily assume that other assumptions 
are fulfilled: age effects may be strong and star-forming regions 
may be individually resolved (Kennicutt & Evans 2012; Krui- 
jssen & Longmore 2014). This means that care must be taken 
when interpreting the SFR. In this case it may be preferable to 
interpret the SFR as a proxy for the local radiation field intensity. 
The dust emission may come from heating by local old stellar 
populations or because of heating by energic photons emitted 
by stars in neighbouring pixels rather than being driven by local 
massive stars. 

3.2. Comparison between monochromatic and hybrid SFR 
estimators 

We now compare popular monochromatic SFR estimators in the 
FUV, Ha (both uncorrected for the attenuation), 24 pm, 70 pm, 
and 100 pm bands with SFR(HQ'-r24 pm), which we take as 
the refererence, to understand how their relation changes with 
the scale considered. We have selected SFR(HQf-r24 pm) over 
SFR(FUV-r24 pm) as we will see in Sect. 6, at local scales the 
24 pm and the Ha are more closely related. We should note 
that in this study we are not so much interested in the absolute 
SFR, which we cannot compute reliably at all scale, as in the 
consistency of SFR estimators with one another and their rela- 
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live variations with spatial scale. These relative variations bring 
us important information on star formation in M33. In addition, 
if different estimators give systematically different results, this 
shows that they cannot all be simultaneously reliable. The re¬ 
lations between the various aforementioned SFR estimators are 
shown in Fig. 3. 

We first observe that monochromatic SFR estimates are well 
correlated with the reference estimates (0.67 < p < 1.00, with p 
the Spearman correlation coefficient). There is a rapid increase 
of the correlation coefficient up to a scale of 150-200 pc for all 
estimators but Ha. Beyond 200 pc, IR estimators show a regu¬ 
lar increase. The FUV correlation coefficient remains relatively 
stable until a scale of 1700 pc and then rapidly increases. The 
Ha estimator globally shows little variation with scale. At scales 
beyond 2 kpc, all estimators are strongly correlated with the ref¬ 
erence one. 

However, if they are all well correlated, this does not nec¬ 
essarily mean that they provide consistent results. In the right 
panel of Fig. 3, we show the mean offset and the dispersion 
between monochromatic SFR estimators and the reference one. 
The FUV, Ha, and 24 pm estimates are lower than the refer¬ 
ence one. This is naturally expected for Ha because it is part 
of the reference SFR estimator. The FUV being subject to the 
attenuation will also naturally yield lower estimates. The am¬ 
plitude of the offset at 24 pm (0.14 dex at 33 pc to 0.10 dex at 
2084 pc) can be more surprising as the 24 pm estimator used 
here is non-linear to take into account that only a fraction of 
photons are attenuated by dust. This is probably due to a metal- 
licity effect. Magrini et al. (2009) measured the metallicity of 
M33 Hii regions at 12 -r \ogO/H = 8.3, placing it near the 
limit between the high (12 -r log OjH > 8.35) and intermediate 
(8.00 < \2 + logOjH < 8.35) metallicity samples of Calzetti 
et al. (2007). In turn, intermediate metallicity galaxies show 
some deficiency in their 24 pm emission relative to higher metal¬ 
licity galaxies. This is due to reduced dust content of the ISM 
which increases its transparency (Calzetti et al. 2007). 

If we compare the SFR at 24 pm and 70 pm, the relative 
offset ranges from 0.57 dex at 33 pc to 0.53 dex at 2084 pc. 
Such a discrepancy has several possible origins. First, these in¬ 
frared estimators have been determined only for a limited range 
in terms of SSFR. Li et al. (2013) computed their estimators 
for -1.5 < logESFR < 0.5 M© yr“^ kpc“^. The 24 pm estima¬ 
tor of Calzetti et al. (2007) benefited from a much broader range: 
-3.0 < logESFR < 1.0 M© yr“^ kpc“^. If we consider only the 
definition range of SFR(70 pm), at the finest pixel size, the dis¬ 
agreement between SFR(70 pm) and SFR(H(r-r24 pm) is not as 
strong. Another possible source of disagreement lies in the scale 
at which estimators have been derived. Indeed, increasing the 
pixel sizes means averaging over larger regions and including a 
larger fraction of diffuse emission. Li et al. (2013) determined 
their estimators on two galaxies at a scale of about 200 pc. They 
estimated that 50% of the emission at this scale comes from dust 
heated by stellar populations unrelated to the latest episode of 
star formation. Yet even if this diffuse emission were exclusive 
to the 70 pm band, this would not be sufficient to explain the 
full extent of the offset. Calzetti et al. (2007) combined data of a 
much more diverse sample of galaxies at a physical scale ranging 
from 30 pc to 1.26 kpc, averaging out specificities of individual 
galaxies. To gain further insight on these differences, we will 
examine in detail the origin of dust emission at different scales 
in Sect. 4. 

As a concluding remark, these discrepancies must serve as a 
warning when using SFR estimators. Their application beyond 
their validity range in terms of surface brightness, physical scale. 


metallicity may yield important biases. This is especially impor¬ 
tant when applying SFR estimators on higher redshift galaxies 
as their physical properties may be more poorly known. 

4. Understanding dust emission to measure the 
SFR at different scaies 

4.1. What the infrared emission traces from 24 jim to 100 jim 

To understand what the emission of the dust traces at which scale 
and under which conditions, we examine the change in the rela¬ 
tive emission at 24 pm, 70 pm, and 100 pm. To ease the compar¬ 
ison, we first convert the luminosity surface densities into SFR 
using the linear estimators of Rieke et al. (2009) at 24 pm and Li 
et al. (2013) at 70 pm and 100 pm. We stress that we are not in¬ 
terested here in the absolute values of SFR but only their relative 
variations. These linear estimators only serve to put luminosities 
on a comparable scale. 

In Fig. 4 we compare the dust emission at 24 pm, 70 pm, and 
100 pm from 33 pc to 2084 pc. In general there is an excellent 
correlation between the emission in these three bands across all 
scales (0.90 < p < 0.98). Unsurprisingly the luminosity of indi¬ 
vidual regions in all bands also varies with ESFR. When examin¬ 
ing relations at a scale of 33 pc, we find that there is a systematic 
sub-linear trend between shorter and longer wavelength bands. 
For higher luminosity surface densities, Ly(24 pm) is stronger 
relatively to Ly(70 pm) than what can be seen at lower Ly(24 pm) 
or Ly(70 pm). The same behaviour is clearly observed when 
comparing Ly(70 pm) with Ly(100 pm). Interestingly, when go¬ 
ing towards coarser resolutions this trend progressively disap¬ 
pears, and at 2084 pc the relations between the various bands 
appear more linear. The important aspect to note is not so much 
that the dispersion diminishes with coarser spatial scales but that 
there is a progressive transition from a non-linear relation to a 
linear relation. This phenomenon could be due to the progressive 
mixing of diffuse and star-forming regions. 

To understand how the relative infrared emission varies with 
the spatial scale, we compare the observed dust at 24 pm, 70 pm, 
and 100 pm with the model of Draine & Li (2007). We refer to 
Rosolowsky et al. (in prep.) for a full description of the dust 
SED modelling of M33 with the Draine & Li (2007) models. In 
a nutshell, the emission of the dust is modelled by combining 2 
components. The first component is illuminated by a starlight 
intensity Umin, corresponding to the diffuse emission. The other 
component corresponds to dust in star forming regions, illumi¬ 
nated with a starlight intensity ranging from Umin to Umax fol¬ 
lowing a power law. We consider all available values for Umin, 
from 0.10 to 25. Following Draine et al. (2007), we adopt a fixed 
Umax = 10^. The fraction of the dust mass linked to star-forming 
regions is y, and as a consequence 1 - y is the mass fraction of 
the diffuse component. We consider y ranging from 0.00 to 0.20 
by steps of 0.01. Because M33 has a sub-solar metallicity, we 
adopt the so-called MW3.1_30 dust composition, which corre¬ 
sponds to a Milky Way dust mix with a PAH mass fraction rel¬ 
ative to the total dust mass of 2.50%, lower than the Milky Way 
mass fraction of 4.58%. We compare this grid of physical mod¬ 
els to the observations in Fig. 5 for a resolution of 33 pc, and at 
a resolution of 260 pc. 

We see that the parameters space spanned by the grid of 
models reproduces the observations very well except for a frac¬ 
tion of points at low 24-to-70 and 70-to-100 ratios, for which 
even models with y = 0 fail. Most points are concentrated in 
regions with simultaneously low values for y and Umin, which 
also correspond to low SFR estimates. Regions at higher SFR 
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Fig. 3. Comparison of monochromatic SFR estimators with the reference SFR(Hcr+24 jjim) estimator versus the spatial scale. The colour 
indicates the monochromatic band: FUV (blue), Her (cyan), 24 {im (green), 70 fim (magenta), and 100 fim (red). The left panel shows the 
correlation coefficient between the monochromatic and reference estimators. The right panel shows the mean offset (solid line) and the standard 
deviation (shaded area) of the difference between the estimators. 
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Fig. 4. Relations at 33 pc (left), 260 pc (centre), and 2084 pc (right) of Lv(70 jjim) versus Lv(24 jjim) (top), and Lv(100 jjim) versus Ly(70 yum) 
(bottom). All luminosities have been multiplied by a constant factor corresponding to a linear SFR estimator (a 24 = 2.04 x 10"^^ M© yr"^ 
ayo = 5.89 X 10"^^ M© yr“^ and aioo = 5.17 x 10"^^ M© yr"^ in order to put them on a similar scale. The colour of each point indicates 
ESFR following the colour bar at the right of each row. The red line indicates a one-to-one relation. We see the non-linear relations between the 
luminosities in different bands. These non-linearities are particularly apparent at the finest pixel scales. At coarser scales, the relations appear 
more linear, which is probably due to a mixing between diffuse and star-forming regions. 
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Fig. 5. 70-to-100 versus 24-to-70 flux density ratios for each pixel at a resolution of 33 pc (left) and 260 pc (right). The colour of each symbol 
corresponds to ESFR, according to the colour bar on the right. The grid represents the Draine & Li (2007) models, with the MW3.1_30 dust 
composition, 0.10 < < 25, and 0.00 <y < 0.20. The red dashed lines indicate the locus corresponding to the one-to-one relations shown in 

Fig. 4. The 3-cr uncertainties are shown in the bottom right comer. We see that the 24-to-70 ratio is well correlated with y and ESFR, especially 
at 33 pc. At 260 pc, due to mixing between diffuse and star-forming regions, excursions in y are strongly reduced. Note that when considering 
the galaxy as a whole, a large fraction of the emission is due to the handful of luminous regions rather than the larger number of faint regions. 


seem to have a higher value for and there is a clear trend 
with 7 , strongly star-forming regions having a larger y. In other 
words, this means that the relative increase of the 24 jjm emis¬ 
sion compared to the 70 yum one that we saw in Fig. 4 is likely 
due to the transition between a regime entirely driven by the dif¬ 
fuse emission and a nearly complete lack of dust heated in star¬ 
forming regions (0.00 < y < 0.01), to a regime with a strong 
contribution from dust heated in star-forming regions. When 
the resolution is coarser the emission from star-forming regions 
is increasingly mixed with the emission from dust illuminated 
by the diffuse radiation field, reducing the excursions to large 
values of y required to have a strong emission at 24 yum com¬ 
pared to the emission at 70 yum. If we assume that on average 
in star-forming galaxies y = \ - 2% (e.g. Draine et al. 2007), 
a significant fraction of the luminosity at 70 yum comes from 
star-forming regions. Considering a resolution of 33 pc, these 
values of y correspond typically to regions with log ESFR > -2 
to -1.5 M© yr“^ kpc“^. The 70 yum luminosity contributed by 
regions brighter than log ESFR = -2 and -1.5 is 58% and 27% 
respectively. This is consistent with what we would expect from 
Fig. 5 as a small fraction of pixels with a high ESFR contributes 
a large fraction to the total luminosity compared to the more nu¬ 
merous but much fainter pixels. 

We can also understand the observed trends by examining the 
physical origin of dust emission in relation to the SFR. At high 
SFR, the emission at 24 yum and 70 yum is caused by dust at the 
equilibrium and by a stochastically heated component. In low 
SFR regions only the stochastically heated component remains 
at 24 yum, contrary to what occurs at 70 yum (see in particular Fig. 
15 in Draine & Li 2007). This means that the 24 yum emission 
should drop more quickly than the 70 yum emission with decreas¬ 
ing SFR. This accounts for the difference in behaviour seen in 
Fig. 5. The preceding explanation for M33 seems consistent with 
the findings of Calzetti et al. (2007, 2010) who have studied this 
problem in great detail. Combining several samples totalling al¬ 
most 200 star-forming galaxies, Calzetti et al. (2010) also found 
a clear positive correlation between the measured SFR and the 
24-to-70 ratio. 


4.2. Impact of the scale on the measure of the SFR from 
monochromatic infrared bands 

4.2.1. Computation of SFR scaling relations 

The determination of the SFR is paramount to understanding 
galaxy formation and evolution. Initially, such estimates in the 
mid- and far-infrared were limited to entire galaxies due to the 
coarse resolution of the first generations of space-based IR in¬ 
struments. Spitzer has enabled the computation of dust emission 
in galaxies at a local scale in nearby galaxies (e.g., Boquien et al. 
2010a). Thanks to its outstanding resolution, Herschel has en¬ 
abled such studies at the peak of the emission of the dust at ever 
smaller spatial scales (Boquien et al. 2010b, 2011; Galametz 
et al. 2013). But such a broad and homogeneous spectral sam¬ 
pling represents an ideal case. More commonly, just one or a 
handful of infrared bands are available at sufficient spatial reso¬ 
lution. It is therefore important not only to be able to estimate 
the SFR from just one or few IR bands but also to understand 
how this is dependent on the spatial scale. 

To do so, we simply determine at each resolution the scaling 
factor Cband between a given band and ESFR from the combina¬ 
tion of Hof and 24 pm. This is done by carrying out an orthogonal 
distance regression using the odr module from the scipy python 
library on the following relation: 

log ESFR = log Cband + log Sband • (1) 

In order to examine the difference between intense and quiescent 
regions, at each resolution we have also separated the regions 
into 4 bins in addition to fitting the complete sample: top and 
bottom 50%, and top and bottom 15%, in terms of ESFR from 
the combination of Ha and 24 pm. The most extreme bins en¬ 
sure that we only select the most star-forming (top) or the most 
diffuse (bottom) regions in the galaxy. 

4.2.2. Dependence of SFR scaling relations on the pixel size 

The dependence of the scaling factors on resolution at 8 pm, 
24 pm, 70 pm, and 100 pm is presented in Fig. 6. 


Article number, page 8 of 15 





































M. Boquien et aL: Measuring star formation with resolved observations: the test case of M33 


I 

© 

hJ 


© 

2 


p-i 

CD 


o 





1.8 

1.6 

1.4 

1.2 

1.0 

0.8 

0.6 

0.4 


le-9 


4.0 


3.5 


3.0 


2.5 


1.0 


Sieke_eLaL C2(K)QJ 


le-9 


Spatial scale [pc] 


Spatial scale [pc] 


I 

© 


I 

Sh 


© 


2 

S 

p-i 

CD 


O 




500 1000 1500 2000 

Spatial scale [pc] 


Fig. 6. Scaling coefficients from the luminosity in infrared bands to ESFR versus the pixel size, at 8 jum, 24 jum, 70 jum, and 100 jum, from the top 
left corner to the bottom right corner. The blue line indicates the value of the scaling factor when taking into account all pixels detected at a 3-cr 
level in all six bands. The red (respectively green) line indicates the scaling factor when considering only pixels with a ESFR higher (resp. lower) 
than the median ESFR at a given resolution. The cyan and magenta lines represent regions in the top and bottom 15% in terms of ESFR. The 
shaded areas of the corresponding colours indicate the 1-cr uncertainties. The horizontal dashed line at 24 jum (resp. 70 jam) indicates the scaling 
factor determined by Rieke et al. (2009) (resp. Calzetti et al. 2010) for entire galaxies. The crosses for the 70 jum and 100 jum bands indicate the 
scaling factor determined for individual galaxies at a scale of 200 pc (Li et al. 2013) and 700 pc (Li et al. 2010). The squares indicate mean values 
over several galaxies. The empty squares denote that no background subtraction was performed. 


Description of the scaling relations It clearly appears that re¬ 
gions with strong and weak ESFR have markedly different scal¬ 
ing factors and a different evolution with pixel size. Compared 
to the entire sample, at 33 pc the scaling factor for the 50% (re¬ 
spectively 15%) brightest pixels is higher by a factor 1.06 to 1.16 
(resp. 1.12 to 1.55). Conversely, the scaling factor for the 50% 
(resp. 15%) faintest pixels is lower by a factor 0.71 to 0.84 (resp. 
0.56 to 0.74). When increasing the pixel size from 33 pc to 
2084 pc, the scaling factor for pixels with a weak ESFR strongly 
increases. On the opposite, the scaling factor for pixels with a 
strong SSFR generally shows a slightly decreasing trend. From 
a typical scale of 400 pc to 1200 pc, depending on the infrared 
band, there is no significant difference in the scaling factors be¬ 
tween pixels with weak and strong SSFR. 


Impact of the relative fraction of diffuse emission As we have 
already explained, our reference SFR estimator combining Ha 
and 24 jim is unfortunately not perfect because it is also sensi¬ 
tive to diffuse emission that may or may not actually be related 
to star formation. We now consider only the 15% brightest pix¬ 


els at 33 pc. They most likely correspond to pure star forming 
regions with little or no diffuse emission. Conversely, the 15% 
faintest pixels will be almost exclusively made of diffuse emis¬ 
sion with little or no local star formation. That way the scaling 
factor will be higher for the former compared to the latter. If we 
move to coarser resolutions, individual pixels will increasingly 
be made of a mix of star-forming and diffuse regions such that 
the brightest and faintest regions will be less different at 2084 pc 
than they are at 33 pc. This naturally yields increasingly similar 
scaling factors that progressively lose their dependence on the 
intensity of star formation. In other words, this means that on a 
scale larger than roughly 1 kpc, monochromatic IR bands from 
8 Jim to 100 jim may be as reliable for estimating the SFR as the 
combination of Ha and 24 jim. This scale is probably indicative 
of the typical scale from which there is always a similar fraction 
of diffuse and star-forming regions in each pixel, bright or faint. 
This scale is likely to vary depending on intensity of star forma¬ 
tion in a given galaxy and on its physical propeties. This aspect 
should be explored in a broader sample of spiral galaxies. We 
also have to mention that this result is also affected by the trans¬ 
parency of the ISM as we will see below, or by non-linearities 
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that are not accounted for here. For instance, in intense star¬ 
forming regions the 8 yum emission may get depressed because 
of the PAH destruction by the strong radiation field (Boselli et al. 
2004; Helou et al. 2004; Bendo et al. 2006), or because the 8 /im 
has a strong stochastic component, proportionally more impor¬ 
tant than at 24 yum. These processes can induce a non proportion¬ 
ality between the 8 yum emission and the SFR. Finally, we note 
that the difference in the scaling factor between the faintest and 
the brightest bins is minimal at 24 yum. This is most likely due to 
the fact that the 24 yum emission affects both sides of Eq. 1 . 

Comparison with the literature When comparing the scaling 
factors determined in M3 3 with those determined in the literature 
from both individual star-forming regions in galaxies and entire 
galaxies, we find instructive discrepancies. At a scale of 200 pc, 
the scaling factors at 70 yum determined by Li et al. (2013) for 
NGC 5055 and NGC 6946 are systematically higher. As dis¬ 
cussed in the aforementioned article, this may be due to back¬ 
ground subtraction. Indeed their study is based on the selection 
of individual H ii regions, allowing for the subtraction of the lo¬ 
cal background, which is not easily doable with accuracy when 
carrying out a systematic pixel-by-pixel analysis like we are do¬ 
ing in this article. Without background subtraction, they obtain 
a scaling factor that is very similar to the one we find when se¬ 
lecting pixels with a strong SFR. A similar study carried out at a 
scale of 700 pc by Li et al. (2010) leads to a similar result. 

When we compare our scaling factors to the ones obtained 
on entire galaxies at 24 yum by Rieke et al. (2009) and at 70 yum 
by Calzetti et al. (2010) there is a clear discrepancy, their scal¬ 
ing factors being lower. Because we see little trend with pixel 
size at larger scales, it appears unlikely that the scaling factor 
will diminish strongly at scales larger than 2084 pc. A possible 
explanation is that this could be due to the increased ISM trans¬ 
parency in M33. In other words, this could be because a smaller 
fraction of the energetic radiation emitted by young stars is re¬ 
processed by dust into the infrared. In the case of M33, about 
75% of star formation is seen in Ha and only 25% in the infrared. 
Indeed, Li et al. (2010) found a trend of the scaling factor with 
the oxygen abundance, with more metal-poor galaxies having a 
higher coefficient. If we consider the relation Li et al. (2010) 
find between the oxygen abundance and the scaling factor, the 
change in the coefficient from 12 -l- log 01H ^ 8.3 (for M33) 
to 12 -h log OjH ^ 8.7 (for the sample of Calzetti et al. 2010), 
would explain the observed discrepancy. At the same time, we 
notice that the discrepancy with Rieke et al. (2009) at 24 yum 
is more important than with Calzetti et al. (2010) at 70 yum. 
This is expected because the former sample is made of the most 
deeply dust-embedded galaxies ([ultra] luminous infrared galax¬ 
ies), contrary to the latter one which is made of galaxies that 
are more transparent at short wavelength. We should however 
note that at a given metallicity, Li et al. (2010) find an impor¬ 
tant dispersion. This is exemplified by the case of NGC 5055 
and NGC 6946, which despite having very similar metallicities 
yield very different scaling factors. In addition, a galaxy like 
the Large Magellanic Cloud which has a metallicity similar to 
that of M3 3 has a scaling factor similar to that of galaxies with 
12 -r log OjH ^ 8.7, perhaps because it has been calibrated with 
HII regions, with diffuse emission having been subtracted, but 
accounting only for the obscured part of star formation (Lawton 
et al. 2010; Li et al. 2010). A dedicated study to disentangle the 
respective impact of the metallicity and the diffuse emission on 
the scaling factors at various scales would be required to fully 
understand this point. 


5. Obscured versus unobscured star formation 

5.1. FUV and Ha attenuation in M33 

Because of the dust we only see a fraction of star formation in 
the UV or Ha. Lollowing Kennicutt et al. (2009), hybrid SLR 
estimators allow us to easily compute a proxy (noted for the 
attenuation (noted A) of the UV and Ha fluxes. 

^Fuv = 2.5 log [SLR(LUV -r 24 yum)/SLR(LUV)], (2) 

Jiua = 2.51og[SLR(Ha-h24yum)/SLR(Hof)]. (3) 

We can also write this more directly in terms of luminosities: 

^Fuv = 2.51og[l+kFuv-24XL(24yum)/L(LUV)], (4) 

Jliia = 2.51og[l+kHa -24 xL(24jum)/L(Ha')], (5) 

with kbandi-band 2 defined as in Sect. 3.1 and Table 1. These ex¬ 
pressions can also be written equivalently in terms of surface 
brightnesses. Before going further, we should keep in mind that 
these estimators have been defined for star-forming regions and 
may not provide us with accurate estimates outside of their defi¬ 
nition range. 

Because the attenuation increases with decreasing wave¬ 
length, the attenuation in the LUV is higher than in the opti¬ 
cal. Lor instance, if we consider the Milky Way extinction curve 
of Cardelli et al. (1989) with the O’Donnell (1994) update, for 
Ay = 1, Ahqt - 0.8 and Afuv - 2.6. However, nebular emis¬ 
sion is more closely linked to the most recent star formation 
episode, and therefore to dust, than the underlying stellar con¬ 
tinuum. As a consequence, the Ha line is actually more attenu¬ 
ated compared to the stellar continuum at the same wavelength 
than what we could expect from the extinction by a simple dust 
screen affecting both components the same way (Calzetti et al. 
1994, 2000; Chariot & Lall 2000). In reality this differential at¬ 
tenuation strongly depends on the geometry between the dust 
and the stars as well as on the star formation history. Given the 
broad range of physical conditions and scales in M3 3, we can ex¬ 
pect the attenuation law between Ha and the LUV band to vary 
strongly across the galaxy and across scales. Such variations 
would provide us with useful information on the effective atten¬ 
uation curve between those two popular star formation tracers. 
The relation between Ha and LUV attenuations as a function of 
ESLR and the specific SLR (sSLR, the SLR per unit stellar mass) 
is shown in Lig. 7. 

We find that on average, the attenuation in M33 is relatively 
low for a spiral galaxy. There are peaks of attenuation reaching 
2.5 mag in the LUV band at a resolution of 33 pc, but when we 
consider large sections of the galaxy at 2 kpc scales the typical 
attenuation is around 0.6 mag in the LUV band and 0.4 mag in 
Ha, making M3 3 mostly transparent in star formation tracing 
bands on large scales. While this is lower than the typical LUV 
attenuation in nearby spiral galaxies (Boquien et al. 2012, 2013), 
it is consistent with previous findings in M33 (Tabatabaei et al. 
2007; Verley et al. 2009). This difference compared to local spi¬ 
rals is probably due to the more metal-poor nature of M3 3. 

Overall, we find that at the finest resolution, regions in M3 3 
span a broad range in terms of absolute and relative attenua¬ 
tions in LUV and Ha. This does not appear to be due to ran¬ 
dom noise though as the locus of the regions appears structured 
according the the intensity of star formation. Regions with in¬ 
tense star formation as traced by the combination of Ha with 
24 yum, tend to have a high JTIfuv compared to Jiua- This is es¬ 
pecially visible at the finest spatial resolution. Intense star form¬ 
ing regions such as NGC 604 show a peak in ^fuv whereas 
no particular increase is seen in Jiua- If we select all pixels 
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Fig. 7. Attenuation in the FUV band versus the attenuation in Her. The colour of each point indicates ESFR (upper row) or the sSFR (lower row), 
following the bar to the right. In the bottom row, the number of regions N and the Spearman correlation coefficient p are indicated. To compute 
the sSFR, the stellar mass in each region was computed from the 3.6 pm emission using the linear conversion factor of Zhu et al. (2010). The 
red line shows the one-to-one relation. The black, magenta, and cyan lines represent the attenuation for a starburst, a Milky Way, and an LMC 
average curve with differential reddening (/ = E{B - V)continuuml E{B - V)gas = 0.44, solid, with E(B - V)continuum being the reddening between 
the V and B bands of the stellar continuum and E{B - V)gas being that of the ionised gas) and without (/ = 1, dashed). For the starburst relation, 
we assumed that even though the stellar continuum follows the starburst curve, the gas still follows a Milky Way curve. Note that the black and 
cyan solid lines are nearly on top of one another. We find that at the finest resolution, there is a broad range in terms of differential reddening. 
Intense star-forming regions have little differential reddening, whereas diffuse regions on the contrary present a strong differential reddening. At 
coarser scales, the averaging between diffuse and star-forming regions yields a differential reddening that is similar to that of starburst galaxies. 
The overall shape of the attenuation law is however only weakly constrained and may vary across the galaxy. 


with ESFR(HQf + 24 pm) >0.1 M© yr ^ kpc ^ at 33 pc, we find 
(^FUv/^Hor) = 3.94 + 1.45, versus (J^fuv/^Hqt) = 1-81 ± 1.11 
for less active regions. As the resolution becomes coarser, ex¬ 
cursions in attenuation become more moderate and the range 
covered in terms of FUV and Her attenuations becomes much 
smaller. At the coarsest resolution, ^fuv and show little 
scatter and they are consistent with a starburst or a Milky Way 
law with a differential reddening (see Sect. 5.2) between the stel¬ 
lar continuum and the gas. What probably happens is that at 
coarser resolutions, intensely star-forming regions and quiescent 
regions merge together, decreasing the dynamic range in terms 
of attenuation properties. At the coarsest resolution, all regions 
have broadly similar properties, which is why they all have sim¬ 
ilar attenuation laws. We detail this aspect in Sect. 5.2. 

Finally, we should also mention the possibility that there 
is a change in the intrinsic extinction laws because of changes 
in the dust composition. Regions at low ESFR are located in 
the outskirts of the galaxy. However this is probably a mi¬ 


nor effect. M3 3 has a very modest metallicity gradient of 
-0.027 + 0.012 dex/kpc (Rosolowsky & Simon 2008). As we 
can see in Fig. 7, a variation of the differential reddening has 
a much stronger effect than a change in the intrinsic extinction 
curve from the Milky Way to the LMC average. 

5.2. Variations of attenuation laws with scale 

At first sight, these variations may seem at odds with the now 
well established picture of differential attenuation between the 
gas and the stars in galaxies (Calzetti et al. 1994, 2000; Char¬ 
iot & Fall 2000). However, this description was conceived in 
the particular context of starburst galaxies and may not apply 
directly to resolved and more quiescent galaxies. Let’s first 
consider M33 at a resolution of 33 pc. As mentioned earlier, 
a low value for SSFR actually corresponds to diffuse emission 
with at most very little local star formation. In this environ¬ 
ment, because gas is intimately linked with dust, the Ha radi- 
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ation always undergoes some attenuation. However, the stellar 
emission may be relatively attenuation-free as it is not partic¬ 
ularly linked to dust, depending on the actual geometry. This 
would explain the relatively shallow effective FUV-Hof attenu¬ 
ation curve that is normally seen in starburst galaxies. Now, if 
we consider star forming regions, the FUV-emitting stars will on 
average be younger and still closely linked to their birth cloud, 
hence undergoing a much higher attenuation than in diffuse re¬ 
gions. Because Ha is always linked to dust, the increase of the 
attenuation is not as strong. If we now consider coarser resolu¬ 
tions, we increasingly mix diffuse and star-forming regions. At 
a local scale J^fuv is on average much larger than Jiua in star¬ 
forming regions but more comparable in diffuse regions, as we 
have seen above. This means that at a global scale, the effective 
Ha-UV attenuation curve should be shallower than intrinsic ex¬ 
tinction curves. This agrees with what we see at a scale of 2 kpc, 
(^FUv/^Ha) = 1-46 + 0.24. 

Using FUV to FIR broadband data on a sample of nearby, 
resolved galaxies at a typical scale of 1 kpc, Boquien et al. 
(2012) found hints of an evolution of the attenuation curve of 
the stellar continuum, with the age of star-forming regions, from 
a starburst-like curve in young regions to LMC-like curves in 
older regions. A consistent result was found on the scale of 
entire galaxies by Kriek & Conroy (2013). They showed that 
0.5 < z < 2.0 galaxies with a high sSFR have a shallower atten¬ 
uation curve. If we assume that a high SSFR is an indication of a 
young age, this would appear to be opposite of the trend we see 
in M3 3. However, a direct comparison is not straightforward be¬ 
cause here we are comparing the nebular attenuation to the stel¬ 
lar continuum attenuation, and with measurements at only two 
wavelengths. In other words, we are looking at the difference 
between the gas and the stellar attenuation curves, measuring 
each at a single wavelength only. 

A major and poorly constrained factor that is im¬ 
portant for this comparison is the differential reddening 
we mentioned earlier, which we can write as: / = 

E(B - V)continuum/E(B - V)gas. This Can also be expressed in 
terms of attenuations. Considering that E(B - V) = Ay/Rv, 

f — Av,continuum/Av,ps ^ Ev,gas/Ev,continuum- As We have Stated 

earlier, in diffuse regions FUV-emitting stars are probably more 
weakly linked to the dust than the ionised gas. As such, in 
diffuse regions / may be much smaller than what it is in pure 
star-forming regions where it should be closer to f = 1. This 
means that in diffuse regions the attenuation of the stellar contin¬ 
uum would be much smaller than the attenuation of the nebular 
emission at a given wavelength. Based on a sample of galaxies 
observed by the SDSS, Wild et al. (2011) found that the opti¬ 
cal depth of nebular emission compared to that of the contin¬ 
uum is significantly higher for galaxies at low sSFR. They at¬ 
tributed this to a variation of the relative weight of diffuse and 
star-forming regions. This means that / is smaller in these more 
quiescent galaxies. Similar results have been obtained by Price 
et al. (2014) based on the 3D-HST survey and by Kashino et al. 
(2013) using ground-based spectra of galaxies at z = 1.6. To ver¬ 
ify these results in M33, in Fig. 7 we have also colour coded the 
relation between ^fuv and Jiua as a function of the sSFR. We 
find a result consistent with that of the aforementioned works. 
Regions with a high sSFR have a high value of / whereas re¬ 
gions with a low sSFR have a low value of /. This way, consid¬ 
ering a variation of /, it is possible that the effective FUV-Ha 
attenuation curve would show a different evolution compared to 
the attenuation curve of the stellar continuum emission. Our re¬ 
sults suggest both a variation of / across the galaxy at a given 
scale from diffuse regions to star-forming regions, and a varia¬ 


tion depending on the scale due to averaging of star-forming and 
diffuse regions that have different values of /. At the finest reso¬ 
lution, a range of / is required to explain the observations across 
the galaxy. However, as we go towards coarser resolutions, the 
observations can be explained with / = 0.44. 

5.3. Limits on the determination of the attenuation 

This discussion relies on the assumption that no systematic bias 
is introduced due to the way we compute the attenuation and 
ESFR. If we consider ESFR from FUV and 24 jum rather than 
from Ha and 24 yum, the trends are not as clear. There is 
a fraction of pixels at 33 pc with very low FUV attenuation 
(0 < JIfuw ^0.1) and moderately high ESFR. This probably 
corresponds to regions with a low level of 24 jum and Ha emis¬ 
sion but with strong FUV. That could be the case for instance in 
a region where recently formed clusters have blown away much 
of the dust and the gas of their parent clouds. Such regions were 
found by Relano et al. (2013), especially in the outskirts of M33. 

A specific bias may affect some diffuse regions. The most 
extreme have JTIfuv < which would require a partic¬ 

ularly strong differential attenuation. A close inspection re¬ 
veals that these regions are also relatively fainter in Ha. The 
relatively higher uncertainties would then propagate into the 
attenuation estimates yielding spuriously low In prac¬ 

tice, they could also be affected by very strong age and radi¬ 
ation transfer effects such as the escape of ionising photons, 
which would reduce the local Ha luminosity, independent of 
the actual attenuation underwent by Ha photons. In that case, 
with little Ha compared to the FUV and for the same amount 
of 24 yum, the selected estimators will then naturally overesti¬ 
mate Ahqt. These regions would in reality not present a dif¬ 
ferential attenuation as extreme as could be inferred from our 
estimates. To ensure that these uncertainties on diffuse re¬ 
gions do not affect our results, we have selected only regions 
with ESFR(HQf -I- 24 yum) > 10“^ M© yr“^ kpc“^, hence remov¬ 
ing purely diffuse regions. We find that we still see the clear gra¬ 
dients described in Fig. 7. This means that if in the most extreme 
regions, the differential attenuation is likely to be overestimated, 
there is still a clear variation of the differential attenuation de¬ 
pending on the sSFR. 

The issues we have presented show the sensitivity of such 
an analysis on the selected SFR estimators and the great cau¬ 
tion that must be used when interpreting such results. A promis¬ 
ing way to reduce such potential problems would be to compute 
the attenuation with a full SED modelling for the stellar con¬ 
tinuum and from the B aimer decrement for the nebular emis¬ 
sion. The increasing availability of spectral maps using integral 
field spectrographs (lES), and large multi-wavelengths surveys 
now makes this possible for nearby galaxies (e.g., Sanchez et al. 
2012; Blanc et al. 2013). Recently, Kreckel et al. (2013) have 
used such IPS data on a sample of 8 nearby galaxies, deriving 
the nebular attenuation from the B aimer decrement and the stel¬ 
lar attenuation from the shape of the continuum between 500 nm 
and 700 nm. Interestingly, opposite to our results and that of 
Wild et al. (2011), they find that in diffuse regions the attenu¬ 
ation of the stars increases compared to that of the gas. In the 
most extreme cases, in the V band the stellar attenuation is 10 
times higher than that of the gas. Conversely, for regions with 
ESFR > 10“^ M© yr“^ kpc“^, they converge to / = 0.47, close to 
what we find at the coarsest resolution. The discrepancy at low 
SSFR may be due to systematics in the way the attenuation is 
computed for the diffuse medium, both for the stars and the gas. 
For similar sized regions at high SSFR, the discrepancy is prob- 
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ably due to the fact we measure the continuum attenuation in the 
FUV whereas Kreckel et al. (2013) measure it in the optical. In 
their case, even in star forming regions the continuum emission 
is generally dominated by older stellar populations, which is not 
necessarily the case in the FUV, inducing a different /. This 
effect is probably prevalent mainly at the smallest scales 

To summarise, great care must be used when correcting star- 
formation tracing bands for the attenuation. We have shown that 
there are clear variations of the effective FUV-Ha attenuation 
curves that depend on the sSFR and ESFR, with regions at higher 
SFR having steeper attenuation curves. This is due to a strong 
variation of the differential reddening between the stars and the 
gas. Intense star-forming regions have little differential redden¬ 
ing (/ ^ 1), contrary to more quiescent regions. Finally, there 
is also a strong variation with the resolution, due to averaging 
regions with different physical properties. At the coarsest reso¬ 
lution, the effective attenuation curve is compatible a differential 
reddening of / = 0.44, which is the value for the starburst curve 
for instance. However, it is not possible to discriminate between 
different laws at fixed differential reddening. 


6. Discussion 

We have found that there is an important variation in the dif¬ 
ferential attenuation in M3 3 with both the spatial scale and the 
sSFR. At the same time, it appears that resolution effects become 
small beyond a scale of 1 kpc. In light of these results, we now 
have a better insight into the relation between UV-emitting stars 
and dust in galaxies (Sect. 6.1). They also allow us to understand 
how the measure of the SFR will be affected by high resolution 
observations with upcoming instruments (Sect. 6.2). 

6.1. Constraints on the relative geometry of stars and dust in 
star-forming galaxies 

The actual geometry between the stars and the dust in galax¬ 
ies is undoubtedly complex, but a simple generalised model has 
emerged for starburst and more quiescent star-forming galaxies 
(Calzetti et al. 1994; Wild et al. 2011; Price et al. 2014). These 
descriptions generally rely on a two-component model frame¬ 
work (e.g.. Chariot & Fall 2000): dense star-forming regions 
and a lower density diffuse medium. We will not come back to 
the general descriptions of galaxies that have been discussed in 
detail in the literature (e.g.. Wild et al. 2011). Our multi-scale 
analysis however sheds light on the distribution between FUV- 
emitting stars and the dust at a local scale. 

In diffuse regions we have found that the differential redden¬ 
ing is large. This shows that the FUV-emitting stars are rela¬ 
tively unassociated with dust. This requires that these stars have 
escaped their birth cocoon or that stellar feedback has induced 
a physical displacement between the young stars one hand and 
the gas and the dust on the other hand. Conversely, the nebu¬ 
lar emission is more strongly attenuated. This means that the 
ionised gas is more associated with dust than the stars. Several 
mechanisms can be invoked. First, this emission may originate 
from gas ionised by nearby massive stars or created by ionising 
radiation that has escaped from more distant star-forming re¬ 
gions. Hoopes & Walterbos (2000) hnd that massive stars in the 
held can account for 40% of the ionisation of the diffuse ionised 
gas in M3 3. An alternative is that it comes from Ha photons 
that have travelled a long distance in the plane of the disk before 
being scattered in the direction of the line of sight. The latter 
possibility is less likely as it would locally boost the Ha lumi¬ 


nosity relative to the 24 jum one, thereby reducing the attenuation 
inferred from Eq. 3. 

Conversely, in star-forming regions there is very little dif¬ 
ferential reddening. This suggests that the UV-emitting stars, 
the dust, and the gas are well mixed and follow similar distri¬ 
butions. The actual geometry drives the transformation of the 
extinction curve, which describes the case when there is a sim¬ 
ple dust screen in front of the sources, into an attenuation curve. 
Constraining the geometry would require additional data to at¬ 
tempt to break the various degeneracies affecting the determina¬ 
tion of the attenuation curve. This is a notoriously difficult task, 
especially since the structure of the ISM is much more complex 
than the simple assumptions that are usually made. 

The progressive convergence towards the canonical differen¬ 
tial reddening of / = 0.44 at larger scales shows the impact of 
the distribution of gas and stars at local scales on the galaxy seen 
at coarser scales. But this also shows the danger of assuming 
similar geometries and attenuation curves across ah scales and 
ah regions in resolved galaxies. Assuming a differential redden¬ 
ing different from what it is in reality can lead to errors of a factor 
several on the determination of the attenuation, and therefore on 
the determination of the SFR. In other words: there is no one 
attenuation law that is valid under ah circumstances. However, 
considering regions of at least 1 kpc across strongly limits res¬ 
olution effects to compute the SFR. This is probably due to the 
broad mixing between star-forming and diffuse regions. We will 
explore in the next section when this scale dependence is most 
likely to have an impact in the era of high resolution observa¬ 
tions. 

We present a simplihed graphical description of the relative 
geometry of stars, gas, and dust in diffuse and in star-forming 
regions in Fig. 8. It is conceptually similar to Fig. 8 in Calzetti 
(2001) but at the same time it shows the fundamental difference 
between normal star-forming galaxies and starburst galaxies. 

6.2. Measuring high redshift star formation in the era of high 
resolution ALMA and the JWST observations 

As we have shown across this article, the determination of the 
SFR or the attenuation is not only luminosity dependent, but it is 
also scale dependent. With the recent commissioning of ALMA 
and the launch of the JWST by the end of the decade, it will fi¬ 
nally be possible to carry out highly resolved observations of star 
formation not only in the nearby universe but also well beyond. 
With such opportunities also come the complexities inherent to 
high resolution studies. To examine in which cases the interpre¬ 
tation of the observations may be affected by the resolution, we 
have plotted in Fig. 9 the physical scale that can be reached in 
the UV and in the IR with the JWST and ALMA. 

Unsurprisingly, the highest resolution will be achieved in the 
rest-frame 200 nm, which will allow us to distinguish 500 pc de¬ 
tails ah the way to z = 10. Unless degraded to lower resolution, 
these images may prove problematic to derive reliably the local 
physical parameters. Conversely, the resolution of rest-frame 
8 fim rapidly degrades with increasing redshift, reaching 9 kpc 
at the maximum redshift of z = 2.5. It would still be consider¬ 
ably useful to carry out resolved studies of low-redshift galaxies, 
a resolution of 1 kpc is already achieved at z = 0.15. 

To probe the peak of dust emission at 100 jum with ALMA, 
a baseline of 1 km appears nearly perfect, with the resolution 
only slightly varying around 1 kpc from z = 3.2 to z = 9.9. 
A baseline of only 150 m would only provide us with a much 
coarser resolution of 4 to 5 kpc. An additional complexity not 
taken into account here would be the loss of uv coverage from the 
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FUV emission from massive stars that have broken out of their birth eloud. 
Only a small fraetion of FUV photons is absorbed by dust. 




Ha emission from gas elouds ionised loeally by massive stars that have broken 
out of their birth eloud and radiation eseaping from star forming regions. ^ 

This ionised gas is elosely assoeiated to dust indueing a ^ 

strong differential reddening respeetively to the stellar eontinuum. g 


Infrared emission from diffuse dust elouds. 


FUV emission from massive stars mixed with gas and dust. 

A small fraetion of FUV photons eseape the star-forming region. 
The eseaping photons eontribute to ionising the ISM. 


Ha emission from ionised gas mixed with dust and massive stars. 
A larger fraetion of Ha photons eseape, following a single 
attenuation law respeetively to FUV. 


Infrared emission from dust mixed with gas and massive stars. 


o 

h-l 


Fig. 8. Simplified description of the relative distribution of stars, gas, and dust at a local scale in diffuse regions (top) and in star-forming regions 
(bottom). The FUV emitting stars are shown in blue whereas older stellar populations are in yellow or orange. The clouds of gas and dust are 
symbolised with grey patches. 


lack of short baselines, which would be especially problematic 
at low redshift. This would require complementary observations 
with the Atacama Compact Array. The addition in the future of 
bands 10 and perhaps 11 will allow the use of shorter baselines 
while extending the window to lower redshift galaxies. Band 11 
would be able to detect 100 yum emission down to z = 1. 

Overall, the synergy between ALMA and the JWST is excel¬ 
lent to probe star formation at a well resolved scale while also 
gaining valuable insight into triggering or feedback. The combi¬ 
nation of these instruments will extend to much higher redshifts 
the spatially resolved multi-wavelength studies that can only be 
done on nearby galaxies currently. 

7. Conclusion 

In order to understand how SFR measurements of galaxies de¬ 
pend on the physical scale, we have carried out an analysis of the 
emission of the local group galaxy M33 from 33 pc to 2084 pc. 
We have found that: 

1. Monochromatic SFR estimators can be strongly discrepant 
compared to a reference }la-\-24 jjm estimator. These dis¬ 
crepancies depend on the scale of the study and on ESFR. 
They may be due to be combined effects of the age, the ge¬ 
ometry, the transparency of the ISM, and the importance of 
diffuse emission. 

2. The scaling factors from individual infrared bands to ESFR 
show an important evolution with the physical size, up to a 
factor 2. Star-forming and diffuse regions show a different 
evolution with the spatial scale. There is however a conver¬ 
gence of scaling factors at large scales. 


3. More generally, such variations with the physical scale and 
the discrepancies of the scaling relations compared to those 
obtained from different samples show that it is especially 
dangerous to apply SFR estimators beyond their validity 
range in terms of surface brightness, physical scale, and 
metallicity. This issue is especially important when apply¬ 
ing SFR estimators on higher redshift galaxies as their phys¬ 
ical properties may be more poorly known. This is why we 
make no attempt at deriving a multi-scale SFR estimator as 
it would be strongly tied with M3 3. That being said, carrying 
out studies at a scale coarser than 1 kpc strongly limits reso¬ 
lution effects. Such resolutions will be routinely achieved at 
high redshift with ALMA and the JWST. 

4. Finally, there is a clear change in the differential reddening 
between the nebular emission and the stellar continuum de¬ 
pending both on the physical scale and on ESFR or the sSFR. 
Star-forming regions have nearly no differential reddening 
whereas diffuse regions have a strong differential redden¬ 
ing. Such a change in the reddening is especially visible 
at the finest spatial resolution. At coarser resolutions, the 
differential reddening converges to values compatible with 
the canonical 0.44 value derived for starburst galaxies by 
Calzetti et al. (2000). These results allow us to obtain new 
insights into the relative geometry between the stars and the 
dust at a local scale in galaxies, from diffuse regions to star¬ 
forming regions. 

The maps presented in this article are available in the form 
of a FITS file at the following address: http: //www. ast. cam. 
ac.uk/~mboquien/m33/. 
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Fig. 9. Spatial resolution versus z for a rest frame wavelength of 
200 nm (yellow) and 8 jjim (red) with the JWST, and at 100 jjim with 
ALMA for baselines of 150 m (blue) and 1 km (green). The Planck Col¬ 
laboration et al. (2014) cosmological parameters are adopted. The solid 
lines show at which redshifts the observations can be carried out. With 
the JWST Observations below 200 nm are not possible under z = 2, 
while observations beyond z = 2.5 are not possible at 8 jjm. Conversely, 
the rest frame emission at 100 fim cannot be observed with ALMA be¬ 
low z = 3.2. There are several gaps at longer wavelengths correspond¬ 
ing to the gaps between different ALMA bands. These bands corre¬ 
spond to the ones available for cycle 2. Band 10 will strongly improve 
the capabilities of ALMA to map the main infrared star formation bands 
at moderate redshifts while band 11 will allow us to probe the peak of 
dust emission down to z = L Finally, the hatched area corresponds to 
a resolution smaller than 1 kpc, where there may be a strong impact on 
the measure of star formation. 
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